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ABSTRACT 

Context. Rich clusters offer a unique laboratory for studying the effects of local environment on the morphological transformation of 
galaxies moving from the blue star-forming cloud to the red passive sequence. Due to the high-density, any environmental process 
should be more pronounced there compared to the field population. 

Aims. Ideally, we would like to reconstruct the evolution of a single, hypothetical representative galaxy as it enters the cluster. For 
that purpose wide-field imaging is crucial to probe a wide range of densities and environments (from the core to the outskirts) and 
isolate this way, the different physical processes which are responsible of the migration from the blue-cloud to the red-sequence. 
Methods. We analyze a sample of 9 massive clusters at 0.4 < z < 0.6 observed with MegaCam in 4 photometric bands (g,r,i,z) from 
the core to a radius of 5 Mpc (~ 4000 galaxies). Galaxy cluster candidates are selected using photometric redshifts computed with 
HyperZ. Morphologies are estimated with galSVM in two broad morphological types (early-type and late-type). We examine the 
morphological composition of the red-sequence and the blue-cloud and study the relations between galaxies and their environment 
through the morphology-density relations (T - E) and the morphology-radius relation (T - R) in a mass limited sample (log(M/M ) > 
9.5). 

Results. We find that the red sequence is already in place at z ~ 0.5 and it is mainly composed of very massive (log(M/M e ) > 11.3) 
early-type galaxies. These massive galaxies seem to be already formed when they enter the cluster, probably in infalling groups, since 
the fraction remains constant with the cluster radius. Their presence in the cluster center could be explained by a segregation effect 
reflecting an early assembly history. Any evolution that takes place in the galaxy cluster population occurs therefore at lower masses 
(10.3 < log(M/M e ) < 11.3). For these galaxies, the evolution, is mainly driven by galaxy-galaxy interactions in the outskirts as 
revealed by the T — £ relation. Finally, the majority of less massive galaxies (9.5 < log(M/M Q ) < 10.3) are late-type galaxies at all 
locations, suggesting that they have not started the morphological transformation yet even if this low mass bin might be affected by 
incompleteness. 

Key words. Galaxies: clusters: general - Galaxies: evolution - Galaxies: high-redshift 



1. Introduction 

Studies of galaxy properties and their dependences on environ- 
ment are key to understanding the mechanisms that are respon- 
sible for galaxy evolution since they constitute a natural test for 
cosmological theories of galaxy formation. Evidence has accu- 
mulated in recent years that environmental processes affect both 
the star formation and morphological characteristics of galaxies, 
the endpoint being the present-d ay morphology-density relation 
dDresslerlll98fi iTreu et al.ll2003l) . 

In particular, two fundamental observational results deal 
with the star formation history and the color-mo rphology dis- 
tribut i on of field galaxies : a) the downsizing (e.g. ICowie et al.l 
Il996t iBundv et al.l I2005I) . i.e. a stronger and earlier decline 
of star formation in massive galaxies compared to a slower 
decrease of the star formation rate of t he present blue/late- 
type galaxies and b) the bimodality (e.g. [Strateva et al.l l200"lt 
iBaldrv et al.ll20041) of the galaxy distribution in the color-mass 
space, i.e. a marked segregation in color between the blue cloud 
(low mass, star-forming late-type galaxies) and the red-sequence 
(massive, red, passive early-type galaxies). 

Therefore, the scenario which seems to emerge is that galax- 
ies move from the blue population to the red-sequence at rates 



and cosmic epochs that depend very strongly on their mass 
and environment. However, the nature and timescale of the 
relevant mechanisms that lead this mass-dependent transition 
are still unclear. We need processes to quench the star for- 
mation above a mass threshold such as gas loss caused by 
AGN feedback and/or supernova feedback (e.g . lDi Matteo et all 
120051; iMenci et alj|2006t iNeistein et al"1l2006h . con sumption of 
all av ailable gas in the wake of violent merging (e.g. lFaber et al.l 
|2007|) . ram-pressure stripping/harassment in clusters that turns 
off efficient cold gas accretion. A combination of these pro- 
cesses are normally required in semi-analytical models to en- 
sure that quenching is efficient, permanent, and mass-dependent 
dDe Lucia et al.ll2006l) . 

From that point of view, rich clusters offer a unique lab- 
oratory for studying the effects of local environment on the 
transition between th e blue cloud and the red sequence (e.g. 
iPoggiantT et al. 2006). Their cores have the highest volume den- 
sity of galaxies in the universe so that any environmental de- 
pendence of galaxy formation or evolution processes should be 
more pronounced there, when contrasted with studies of the 
field population. Ideally, we would like to make a significant 
step forward and reconstruct the complete evolution of a sin- 
gle, hypothetical representative galaxy as it enters the cluster 
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potential and interacts with the ICM. How do star-formation and 
morphology evolve as a galaxy enters the cluster? How does it 
move form the blue-cloud to the red-sequence? To answer these 
questions it becomes necessary to assign timescales to the phe- 
nomena observed and to connect afterwards clusters at differ- 
ent redshifts into a temporal sequence. The various processes 
that can take place can be classified in three broad headings 
(see iTreu et alj|2003l for a nice review): 1) galaxy-ICM interac- 
tions, which refers to the interaction of a cluster galaxy with the 
gaseous component of the cluster, 2) galaxy-cluster gravitational 
interactions, including tidal and related dynamical processes and 
3) smaller scale galaxy-galaxy interactions. 

During the past years, much progress have come from the 
connection of Hubble Space Telescope (HST) imaging with deep 
ground-based spectroscopy. Several HST-based studies have 
demonstrated that the blue galaxies seen at z > 0.3 are star- 
forming galaxies which are absent in pre sent day clusters, (e.g. 
iDressler et~ai]|1994t IStanford et~ai1ll995l) It has been proposed 
that these are recent arrivals in the clus ter, probably observe d 
before the removal of their gas reservoir dDe Lucia et alj|2007h . 

The main problem of HST imaging is however its small field 
of view, which makes it very time consuming to obtain a wide 
view of the cluster from the inner core to the outskirts. That is 
the main reason why most of these previous works focus on the 
inner 0.5 Mpc of several clusters which is lower than the virial 
radius (~ 2Mpc). However probing the outskirts of clusters is 
essential to discriminate between the different processes that can 
induce morphological tr ansformation of a galaxy entering the 
cluster. ITreu et alT ( l2003l) have shown indeed that depending on 
the analyzed zone it is possible to isolate one process from an- 
other. Thanks to a wide survey on a single target, they were able 
to show that the most likely processes responsible for the mild 
gradient in the morphological mix outside the virial radius are 
galaxy-galaxy interactions. This later affirmation has been con- 
firmed in more r ecent works. Using a volume-limited sample, 
iPatel et all (120091) argued that the evolution in the galaxies enter- 
ing a cluster at z ~ 0.8 is mainly driven by the local environment. 
In this last work, the authors also showed an interesting trend: 
the dependence of galaxy evolution in clusters with the stellar 
mass. They indeed identify the same downsizing effect as seen 
in the field dCowie et al.lll996tlBundv et al.ll2005l) . with massive 
galaxies having their star forming rate decline faster. 

At this point, it becomes interesting to have a wide view of 
a statistically representative sample of clusters at high redshift 
with robust stellar mass estimates to probe the physical mech- 
anisms responsible of the mass-dependent transition from the 
blue-cloud to the red-sequence. MegaCam, the wide field im- 
ager installed at the Canada France Hawaii Telescope can probe 
in a single pointing a cluster at z ~ 0.5 from the core to ~ 12 
Mpc radius. Problems come now from angular resolution, since 
determining morph ologies of galaxies from this g round-based 
images is not trivial. iHuertas-Companv et al.l d2008l) describes a 
method to estimate robust morphologies on poorly resolved data. 
In that work, the authors were able to prove that, when applied 
to seeing-limited data, the method can provide morphological 
information in two broad classes, with an accuracy close to that 
obtained with space data. 

In the present paper we use this method combined with ac- 
curate photometric redshift estimates in 9 rich clusters at inter- 
mediate redshifts (0.4 < z < 0.6). We show that we are able 
to isolate galaxy cluster members and to estimate the morpholo- 
gies of ~ 4000 galaxies with a fair degree of reliability. We then 
connect them to the cluster environment in three mass ranges to 



probe the mass-dependent transition from the blue-cloud to the 
red-sequence. 

The paper is organized as follows: in §[2] we introduce the 
observational data. In §[3] we describe our method for selecting 
member galaxies based on photometric redshifts. We combine 
the photometric redshift estimate with the shape of the redshift 
probability function (P(z)) function to establish a robust crite- 
ria for selecting galaxies. The physical properties computed are 
described in §[4] We specially focus on the morphology determi- 
nation and on the accuracy quantification. We then describe and 
discuss our main results in § [5] 

Throughout the paper, we assume that the Hubble constant, 
the matter density, and the cosmological constant are Ho = 70 
km s _1 Mpc -1 , Q.,„ = 0.3 and Qa = 0.7 respectively. 

2. Data description 

The data used in this paper are part of an imaging survey of 
X-ray clusters at intermediate redshift (z ~ 0.5), undertaken to 
study the mass distribution in these clusters and to compare the 
X-ray distribution of the intra-cluster gas and the dark matter dis- 
tribution traced by weak lensing. The sample has been selected 
to be representative of the rich cluster population at intermediate 
redshift and is part of the full and unbiased sample observed with 
XMM-Newton as part as a Large Programme (PI: M. Arnaud). 
As the main goal of the optical follow-up was to study in de- 
tails the mass distribution by weak lensing techniques, only the 
hottest and most massive clusters of the sample were selected. 
This ended with a sub-sample of 10 clusters, 9 of which are 
studied in this paper. The last one, RXJ 1347.5-1 145 was dis- 
carded because it was observed in different conditions at CFHT 
and so would reduce the homogeneity of the study. The 9 result- 
ing clusters were observed with the wide-field imager MegaCam 
installed at CFHT, and in 4 photometric bands (g', r', i',z') and 
with integration times of 1600 sec, 7200 sec, 1200 sec. and 
1800 sec. respectively. We asked for a longer exposture time in 
r' because these data are dedicated to the weak lensing mea- 
surement of the mass distribution in the clusters, while the color 
information is restricted to the study of the cluster population. 
We will therefore use r' band data for the morphological analy- 
sis as well (§ 14. U . The total field of view covered by the mosaic 
is one square degree, with a pixel scale of 0.187"(~ 1-1 Kpc at 
z = 0.5). At the average cluster redshift (z ~ 0.5), this corre- 
sponds to a radius up to ~ 12 Mpc, well above the virial radius 
(~ 1 - 2 Mpc, see table [TJ. Consequently, the outskirts of the 
clusters as well as they global environment can also be probed. 
The basic properties of the clusters are summarized in table Q] 
and more details can be found in Foex et al. (in preparation). 

Data reduction was done in a standard way in the Terapix 
cente£] or using the specific tools developped by Terapix (as- 
trometry, stacking and photometric calibratio n). The photomet- 
ric ca talogs were built using SExtractor dBertin & Arnoutsl 
1 19961) in dual mode for multi-color photometry. As illustra- 
tion, we show in Fig. [1] a color composite image of one of the 
clusters. On average, the image quality on the final r' images 
is 0.7 + 0.1" with some variation from cluster to cluster. The 
completeness magnitude of the photometric catalogues is about 
r' = 25, g' = 25, i' = 24 and z' = 23.5. For the details of 
the photometric survey of these clusters, we refer to Foex et al. 
(2009, in preparation). In this paper, we restrict the analysis to 
objects with 18 < r < 23 lying at a radius lower than 5 Mpc. 
As we will show in section $4] this is equivalent to a mass cut 

1 http://www.terapix.fr 
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at log(M»/M ) > 9.5. The magnitude cut is motivated by the 
reliability of the morphological classification (§ 14. It while the 
distance cut is done to reduce the field contaminations (§|3J. 

3. Cluster member selection 

The most unambiguous way to determine cluster membership 
is by means of accurate spectroscopic redshifts. Unfortunately, 
it is too time consuming to obtain high-sepctroscopic complete- 
ness, even to relatively bright limits. He we adopt a method for 
membership deter mination based on photometric redshift esti- 
mates presented in lPello et al.l(l2009l) . Usually, z p /, or based meth- 
ods for determining cluster membership are based on a simple 
cut in redshift, such that a galaxy is considered to be member if 
IZphot—ZdusteA < Azthresh- The problem of this kind of approach is 
that Zthresh can be as high as 0.3, producing high field contamina- 
tions. As a matter of fact, the accuracy of photometric redshifts 
as a function of redshift strongly depends on the filter set and 
its ability to detect strong spectral features, such as the 4000A 
break. In this case, the contiguous spectral coverage from g to 
z is well suited to derive good-quality photometric redshifts for 
galaxies between 0.2 < z < 1.2, i.e. within the redshift range of 
interest for our study. Galaxies below z ~ 0.2 or beyond z ~ 1 .2 
will have poor-quality redshifts and broad redshift probability 
distributions. We therefore expect this population to contami- 
nate the Zphoi distribution in these regions (see Fig.|4|i. To reduce 
this contamination, we combine both the best-fit redshift and the 
information contained in the full redshift probability distribution 
P(z), by computing the integrated probability for the galaxy to 
be in the cluster (P c iust)' 

r-Zdusur+OA 

Pcim = P(z)dz (1) 

Since the contaminating galaxies do not exhibit strong fea- 
tures within the griz spectral domain (i.e. their broad redshift- 
probability distributions avoid the cluster redshift), this selec- 
tion criterion based on the probability distribution should re- 
move most of them. By using a representativ e spectroscopic 
samp le, which is available for MS045 1.5-0305 (lEllingson et al.l 
1 1998b we can calibrate the threshold value Pthresh above which 
a galaxy is considered to be a cluster member. This criteria is 
then extrapolated to the rest of the sample, since all the clusters 
lie in the same redshift range and are observed in similar con- 
ditions (same detector, integration time and set of filters). We 
apply this method to all the objects detected with SExtractor 
within r < 5 Mpc and 18 < r < 23. We use the HyperZ code 
dBolzonella et al.1 [20001) and the 4 available photometric bands 
to compute photometric redshifts (Fig.|4]i. For color indices, we 
use the MAG_APER magnitudes, which give the flux on a fixed 
2.5 arcsec aperture, corrected from zeropoint errors. As there re- 
mains some uncertaintites in the photometric calibration in some 
filters, we add some positive or negative offsets to the zeropoints 
of the different filters determined so that the color-color relations 
of the stars de tected in the images fit those of the Pickles library 
(lPickleslll998h . 

Figure [3] shows the evolution of the rejection of field galaxies 
as well as the completeness of the members selection for differ- 
ent values of Pthresh and for different radii. We clearly see that 
increasing the probability threshold results in an increase of the 
fraction of rejected non-members while the fraction of cluster 
members progressively decreases. Since our main objective is 
to study the effects of the cluster environment on morphology 
and color we want to reduce the number of field contaminations. 




Pthresh 

Fig. 3. Galaxy cluster member selection accuracy for MS 
0451.5-0305 . The different panels show the results for different 
radii. Solid lines indicate the completeness and dashed lines the 
contaminations from non-members as a function of the probabil- 
ity threshold. The vertical dashed line shows the adopted proba- 
bility threshold. 

We therefore consider as galaxy cluster members all the objects 
which have \z p hoi - ZdusteA < 0.1 and P duster > 55%. As shown 
in Fig. [3] this ensures at least 80% completeness and ~ 80% 
rejections up to r < 2Mpc. Above this distance, the complete- 
ness is slightly lower (~ 70%) but the rejection remains high. 
We establish the limit of 5 Mpc as the distance at which there 
is still a reasonable compromise between completeness and field 
contamination, well above the virial radius of the clusters. This 
results in a sample of 3861 galaxies distributed in the 9 clusters. 
We show in figure [2] the spatial distribution of galaxies in the 9 
clusters. 

4. Physical properties 

4.1. Morphology 
4.1.1. Method 

We now discuss the morphological classification of the objects 
belonging to the cluster. Morphologies are determined in the r- 
band images using galSVMq galSVM is a non-parametric N- 
dimensional code based on support vector machines (S VM) that 
uses a training set built from a local visually classified sample. 
dHuertas-Companv et alj|2008l) . Basically we separate galaxies 
in two broad morphological types (late and early type). By late- 
type we mean spiral and irregular galaxies (S/I) and early-type 
galaxies include elliptical and lenticular types (E/S0). 
Broadly, t he employed procedure can be summarized in 4 main 
steps (see lHuertas-Companv et ai]|200 8. for more details): 

1. Build a training set: we select a nearby visually classified 
sample at wavelengths corresponding to the rest-frame of the 

2 http : //www . lesia . obspm . f r/~huertas/galsvm . html 
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Fig. 1. Color image of a 6 x 6 (~2.5Mpcx2.5Mpc) region around the center of MS 0451.5-0305 . The image has been built with 
MegaCam g, r and i filters. 



high redshift sample to be analyzed to minimize morpholog- 
ical k-correction effects. In our case, we want to simulate r- 
band observations at z ~ 0.5. We use therefore an SDSS local 
sample observed in the g band which roughly corresponds to 
the rest-frame wavelength. We then move the sample to the 
proper redshift and image quality, assign a magnitude ran- 
domly selected in the magnitude range [19:23] with a dis- 
tribution representative of the observed luminosity function 
and drop it in a real image to properly simulate the back- 
ground noise. 

2. Measure a set of morphological parameters on the sam- 
ple. In this particular case, the classification is made 



in 7 dimensions, including: ellipticity, asymmetry, gini, 
concentration, M20 and smoothness. For more details on 
how these parameters are exa ctly computed please consult 
iHuertas-Companv et al.1 d2008l) 

3 . Train a support vector based learning machine with a fraction 
of the simulated sample and use the other fraction to test and 
estimate errors. 

4. Classify real data with the trained machine. 
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Table 1. Main properties of clusters in our sample. m lBohringer et all (120001) . [2! lSadat et all d2004l) . [3] Foex et al. 2009, in prep. 
r2oo are computed through a ID weak-lensing analysis, i.e. fit of the radial shear profile and are defined as the radius at which the 
density is 200 times the critical density of the Universe at the cluster redshift. 
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Fig. 4. Photometric redshift distributions for all the clusters. 
Dashed vertical lines show the cluster spectroscopic redshift. 



4.2. Accuracy of classifications: cross-checking with HST 

One of the main challenges of the present work is performing 
a morphological classification on seeing-limited images. Until 
now, all the works involving morphological classifications of 
distant galaxy clusters have been based on HST imaging. We 
need consequently to carefully investigate the accuracy of our 
classification. Fortunately, two of our clusters have published 
data obtained with the space telescope. This is the c ase of MS 
1621.5+2640 dEllingson et alj|1997t ISaintonge et al 
MS 0451.5-0305 dEllingson et al Jl 19981: iMoran et al. 



5005l) 
20071) . 



and 
We 

therefore have a reference to calibrate our morphological classi- 
fication on seeing-limited images. 

MS 0451.5-0305 This cluster has available published results 
on 96 galaxies based on HST i maging which enabl es a second 
calibration of our classification. IMoran et al.l d2007l) performed 
detailed visual morphological classifications using WFPC2 



Fig. 5. Stamps of the same galaxies seen by HST and MegaCam. 
The stamps size is 10 x 10 . For each galaxy, we indicate 
the morphological class obtained from the MegaCam image 
(L=late-type, E=early-type) as well as t he probability of be ing 
early-type and the visual classification of lMorane t al. (2007|) for 
the HST stamps. 



imaging. Figure [6] shows the comparison of our GALSVM-based 
classification with their results. Both classifications are in good 
agreement: a big fraction (~ 90%, 19/21) of early-type objects 
according to HST (i.e. E +S0) fall in the galSVM early-type 
class and 33/41 late-type galaxies (i.e. S+Sc+Sd+I) fall in our 
late-type class. Most of the problems therefore seem to come 
from the intermediate morphology class. Indeed, we are not 
able to separate properly the Sa/b galaxies which seem to ran- 
domly fall in one or the other class. It implies then that we are 
slightly over-estimating the early-type fraction. For illustration, 
we sh ow in figure [5] som e stamps of the same galaxies seen by 
HST (IMoran et alj|2007l) and MegaCam. 
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Fig. 2. Spatial distributions of galaxies in the different clusters. Red circles indicate early-type objects and blue spiral symbols late- 
type galaxies. The symbol sizes code three different mass ranges; big: log(M/M ) > 11.3, intermediate: 10.3 < log(M/M ) < 11.3 
and small: 9.5 < log(M/M G ) < 10.3. Black circles show physical radii: 1, 2 and 5 Mpc. Contour lines indicate density levels: 
D > 2Mpc~ 2 (solid line), D > l.5Mpc~ 2 (dashed line), D > lMpc~ 2 (dashed-dotted line) 



MS 1621.5+2640 ISaintonge etaD (l2005h performed an auto- 
mated morphological classification based on model-fitting of 
258 galaxies. They used the bulge fraction (B/T) as the main pa- 
rameter to split objects into different morphological classes. In 
table [2] we compare the probability of being ea rly-type accord- 
ing to our classification and the bulge fraction of lSaintonge et alj 
d2005l) for the 59 galaxies identified as cluster members and 
within our selected magnitude range: a galaxy with a low bulge 
fraction should have also a low probability of begin elliptical. 
Of course the dispersion of the distribution is considerable but 



we have to keep in mind that we are not measuring the same 
parameters, so points do not necessarily have to lie on the unity 
slope line. Since we are only interested in splitting the galax- 
ies into two rough morphological types, we have to look at the 
fraction of galaxies that are clearly miss-classified, if we con- 
sider the HST-based classification as a reference (which is not 
necessarily true since it is an automated classification as well). 
This includes objects which are clear disc-dominated galaxies 
{B/T < 0.3) but have a probability of being early-type greater 
that 0.5 and vice- versa, galaxies with a bulge fraction bigger than 
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E E/SO SO Sa + b S Sc + d l/M 
HST visua 



Fig. 6. Relat i on be tween the visual morphology computed in 
iMoran et all (120071) and the galSVM type. Filled histograms 
show the fraction of galSVM early-type objects and empty his- 
tograms the fraction of late-type galaxies. 





Pearly > 0.5 


Pearly < 0.5 
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30% ( 11) 


5% (1) 


B/T < 0.3 


10% ( 5 ) 


62% (11) 
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60% ( 25 ) 


33% (6) 



Tabl e 2. Relat i on be tween the bulge fraction computed in 
ISaintonge et alj d2005l) and the probability of being early-type 
computed by galSVM. 



0.7 but low probabilities of being ellipticals. We find that the first 
represent 5 objects out of 41 early-type objects and the second 
one out of 18. We conclude then that our classification is robust 
and in good agreement with the HST one. Note however that we 
might be slightly over-estimating the fraction of elliptical galax- 
ies if compared with HST studies as pointed out in the previous 
comparison. If we indeed consider galaxies with an intermedi- 
ate bulge-fraction (i.e. 0.3 < B/T < 0.7), most of them (25/31) 
fall in the early-type class. It worth noticing here that despite of 
these differences with HST, this work presents 2 important ad- 
vantages: first, the size of the sample (~ 4000 galaxies) makes it 
more robust to miss-classifcations and second, we are using the 
same method to classify galaxies up to 5 Mpc, in all density re- 
gions, which implies that any difference between the center and 
the outskirts must be real. 

4.3. Stellar masses 

Recent works have shown that the relations between galaxies and 
their enviro nment in clusters strongly depend on t he considered 
mass r ange dPatel et alj|2009t IHolden et al.ll2007l) . iHolden et al.l 
d2007l) did not find, for example, any evolution on the early-type 
population with M > 10 10 6 solar masses, using a mass-limited 



sample. They suggested therefore that most of the morpholog- 
ical transformation which occurs from z ~ 0.8 must affect the 
low mass galaxies. It appears therefore interesting to investigate 
the effect of stellar mass on the different analyzed trends in our 
sample. In this work, stellar masses are computed directly from 
t he best-fit t ing m odel of the photometric redshift estimate using 
a lChabrierl (120031) IMF. The mass estimate takes into account the 
age of the model as well as the star forming history. Even though 
we do not have NIR photometry, we consider our mass estimates 
reliable since at z ~ 0.5 the redder filters (i,z) are probing the 
SED around 5000A, (beyond the 4000A break) where the flux is 
mostly coming from evolved stellar populations. 
Mass distributions for all the clusters are shown in figure|7] They 
are clearly bimodal for most of the clusters, with late-type galax- 
ies dominating the low-mass regime (log(M/M G ) < 10.5) and 
the early-type ones dominating the high-mass one log(M/M Q ) > 
11, with of course some cluster-to-cluster variations. We try to 
estimate the mass completeness of our sample using the magni- 
tude cut for the morphological classification (r' < 23), which is 
about 2 magnitudes brighter than the magnitude limit of the cat- 
alogs. If we compare the mass function of the morphologically 
classified sample with that of the full sample up to the magni- 
tude limit (r' = 25), we find that it is complete at the 50% level 
for a mass of log(M/M ) = 9.5 (or M a 2.5 1O 9 M ). In the fol- 
lowing, we will therefore divide the sample in three stellar mass 
bins: 9.5 < log{M/M Q ) < 10.3, 10.3 < log(M/M Q ) < 11.3, 
log(M/M Q ) > 11.3. While the two last bins are not affected by 
incompleteness, the less massive one can be, specially for the 
early-type population which is mostly red and is therefore lost 
earlier in a magnitude limited survey. Nevertheless, we have de- 
cided to include this mass bin since some information can be 
extracted anyway. However, we are aware that this might intro- 
duce some biases when interpreting our results so we will prop- 
erly point out and take into account this point when necessary 
in the following sections. Note also that a rough estimate of the 
apparent magnitude of a L* galaxy at z ~ 0.5 is m* = 20.6 so the 
magnitude limit for the morphological classification corresponds 
to a L* / 10 galaxy. 

4.4. Local density 

The local density has been shown to play a key role in the evolu- 
tion of galaxies, since the efficiency of the different mechanisms 
involved in such transformation often depends on the local den- 
sity. To characterize it for all the galaxies in our sample we com- 
pute the local projected galaxy density, S, at each galaxy loca- 
tion. We use the n-th nearest neighbor definition, with n=10, to 
compute E: 



where d„ is the projected distance of the n-th nearest neighbor. 
The value of n = 10 is commonly used in clusters since it is 
a good measu re of the local ove r-density in this high density 
environments (ICooper et al.ll2005l) . We include in the computa- 
tion, all the cluster members selected with the method presented 
in § [3] To avoid edge biases, we include galaxies a bit further 
of our 5Mpc distance limit. We show in figure [8] the relation- 
ship between radius and local density. We clearly see that local 
density is not a well-defined function of radius, except in the 
central ~ 1 Mpc. We confirm then that by studying the cluster 
outskirts we are able to isolate the effects of over-density from 
the radius , and deduce which one of both effects most influences 
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Fig. 7. Stellar mass histograms for all the galaxies in the 9 clusters. Solid black lines show the histogram for all the galaxies while 
color lines show the histograms separated by morphological type (blue: late-type, red: early-type). 



the morphological evolution of galaxies (§ 15. 31 >. Departures from 
spherical symmetry in the distribution of galaxies in the clusters 
will be enhanced by using the local density instead of the radial 
distance to the cluster center. This is more prominent in clus- 
ter outskirts, where anisotropies of the cluster distribution or the 
presence of in-falling groups are better taken into account. That 
way, we can better evaluate how local effects influence the mor- 
phological evolution of galaxies (§ 15. 3b . 



5. Results and discussion 

5. ) . Composition of the red-sequence 

Mass-color diagrams are a powerful tool to study the star- 
formation (SF) activity of galaxies. Galaxy distribution in this 
plane is bimodal, with passive red galaxies lying in a tight 
mass-color relation called the red-sequence (RS) and blue-star 
forming galaxies presenting a wider range of colors and masses 
and therefore lying in the so-called blue-cloud. This bimodality 
is usually analy sed through the color-magnitude relation (e.g. 
iMei et al1l2009l) . The advantage of the color-stellar mass rela- 
tion is that it displays a less scattered relation, especially for the 
blue cloud. This is can be understood because there is a correla- 
tion between color and M/L ration for late-type galaxies which 
is partly removed in a stellar mass-colo r relation and gives a nar- 
rower blue cloud (e.g. lPatel et al.ll2009l) . In terms of galaxy evo- 
lution in dense environments, as galaxies enter in the cluster, we 
expect that their star-formation processes are quenched by the 




Fig. 8. Local density vs. radius for all the clusters in our sample. 



different interaction effects. Galaxies move progressively from 
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the blue-cloud to the red-sequence, while suffering a morpholog- 
ical transformation from late-types to earlier types. By looking 
at the composition of the red-sequence we can identify which 
galaxies have already suffered the transformation in clusters at 
z ~ 0.5 and then try to isolate the processes responsible of this 
transformation. Nevertheless, this scheme is too simple to inter- 
pret all the diversity of galaxy morphologies and cannot explain 
the status of the "red late-type galaxies" or the "blue early-type" 
ones, but at least part of the bimodality is imprinted by evolution 
processes. 

Figure [9] shows the g-r-stellar mass diagrams in the central 
2 Mpc of all the clusters in our sample. We define whether a 
galaxy belongs or not to the red-sequence based on its deviation 
from the color-mass relation (CMR) shown in figure [9] (black 
solid line). The CMR is derived by performing an iterative fit 
to early-type galaxies with log(M/M&) > 10.5 and rejecting 
the 2o" outliers of the best-fit at each iteration upon conver- 
gence. A galaxy is then said to belong to the red-sequence if 
it lies within the 2cr scatter of the non-rejected galaxies (dashed 
lines). It is worth noticing that despite of the fact that we see 
the same trends in all the clusters, we find that the CMR is 
not universal and that it is subjected to a cluster-to-cluster vari- 
ance over the cluster population with different scatters and abun- 
dances. It is beyond the scope of this paper to analyze in de- 
tails the shapes of the red-sequences but we are planning to 
do it in a coming paper. However, some interesting informa- 
tion can be extracted from these plots. First, they help check- 
ing the mass completeness limit of our sample and verifying that 
red objects (mostly early-type) are indeed lost earlier. As a mat- 
ter of fact, above log(M/M e ) ~ 10.3 we do clearly see the RS 
while below this mass, red galaxies start to become less abun- 
dant. This fact confirms that results obtained in the low mass bin 
(9.5 < log(M/M Q ) < 10.3) might be affected by incompleteness 
and should be analyzed carefully. Second, it is possible to study 
the general composition of the red-sequence. If we first look at 
the morphological distributions, we clearly find out the expected 
trend, i.e. the majority of object lying in the RS are early-type 
objects (76%, table[3]l. This result just confirms that as galaxies 
enter in the cluster and transform into elliptical galaxies, they 
progressively become redder and move to the red-sequence. 
Furthermore, the majority of early-type objects lying the red se- 
quence are massive objects (log(M/M Q ) > 11.3) and as we pro- 
gressively move into lower masses the number of galaxies of 
later morphological types becomes more important. 

These tendancies with the stellar mass suggest that most of 
the evolution in these intermediate redshift clusters is affecting 
low mass galaxies (log(M/M G ) < 11.3), while the most massive 
galaxies seem to be well in place in the red-sequence. The se re- 
sults seem to be in agreement with iDe Lucia et al.l (|2007) who 
also found that, at intermediate redshift, luminous galaxies are 
already in place in the RS and evolve passively, while there is 
a deficit of faint red galaxies, suggesting that these galaxies are 
populating the RS from the blue cloud between z ~ 0.8 and 
today . Furthermore, semi-analytical models (e.g. iMenci et"ai1 
2008) also predict that at a cosmic time of t ~ 2.5 Gyr (z ~ 2.5), 
more than 90% of the stellar mass finally assembled in red- 
sequence galaxies at z ~ 1 is already in place. 

5.2. Radial trends 

The morphologies of galaxies are correla ted with distance from 
the cluster center in nearby clusters (e.g. l Whitmore etaT1ll993h 
and groups (e.g. Postman & Gelled 11984" as well as at highe r 
redshifts (e.g. iAbraham et alJ 119961: Ivan Dokkum et all feOQOl) . 



An increase of the fraction of early-type galaxies towards the 
center is observed in all cases, but the detailed evolution of 
the morphology-radius relation from high redshift remains 
uncertain. Figure [10] shows the fractions of early and late-type 
galaxies at different radii for clusters in our sample. All the 
clusters present similar trends: as detected in previous works, 
we observe a correlation between morphology and distance, 
with elliptical galaxies representing about 80% of the galactic 
populations in the inner parts (d < 0.5 Mpc) and decreasing 
at larger distances. More precisely, we see a steep decrease 
of the early-type fraction in the central Megaparsec which 
goes from ~ 80% to ~ 60% and then a mild decrease in the 
outer parts from ~ 60% to ~ 40% in ag reement with previous 
HST-based works (e.g. iTreu et alJ [20031) . These two different 
regimes are normally explained by different processes. While 
the explanation of the steep increase of the early-type fraction 
in the inner cores requires phenomena related to the cluster po- 
tential, such as tidal interaction or ram-pressure stripping which 
accelerate the morphological transformation, this is not the case 
in the outskirts. As a matter of fact, in the periphery, most of 
the galaxies have never been through the center and therefore 
are free from effects such as ram-pressure stripping. The mild 
decrease in the outskirts might be explain ed by a progressive 
quenching of the star formati on (starvation, Balogh et al.ll2000h 
or after a starburst phase (e.g. lPoggianti et al.lll999l) . 

Another possibility is that the morphology radii relation is 
the consequence of a segregation effect: objects with different 
radii represent galaxies with different early assembly histories 
(through mergers) that keep a cor relation between thei r history 
and the location in the cluster (e.g. lDiaferio et al1l200ll) . 

If this last scenario reveals to be true massive galaxies should 
have had their star formation activity decrease faster than less 
massive galaxies and therefore should be already predominantly 
early-type when they enter the cluster at z ~ 0.5. Figure QT| 
shows the radial trends for different mass ranges, averaged for all 
the clusters in order to have better statistics. The fraction of mas- 
sive galaxies (log(M/M Q ) > 1 1.3) remains almost constant at all 
radii, indeed suggesting that these galaxies were formed before 
they entered the cluster, probably in infalling groups through 
galaxy-galaxy interactions and consequently experience small 
evolution. As a matter of fact, evidence that the source of the 
constituent galaxies of massive clusters ar e groups in the out- 
skirts has also been seen in previous works dDressler et al.f 19971 
iGonzalez et aT1l2Q05h iMoran et al.lf2007h iPatel et al.ll20'09l) . The 
presence of these galaxies in the cluster center can perhaps be 
explained by a segregation effect. 

The general trend of a decrease of early-type galaxies with ra- 
dius has to come therefore from lower mass galaxies, especially 
the intermediate mass range (10.3 < log(M/M e ) > 11.3). Most 
of the evolution in the cluster galaxy population therefore seems 
to be occurring in this stellar mass range. On the contrary, for the 
lowest mass bin in our sample (9.5 < log(M/M Q ) < 10.3), the 
galaxy population is dominated by late-type galaxies as if they 
did not have time to evolve yet. Note that we do not consider 
the very central part of the cluster (d < 0.5 Mpc) because there 
are too few low mass galaxies and numbers are not significant 
enough. We have to be however very cautious when analyzing 
this low mass bin, since as stated in 35.11 it might be strongly af- 
fected by incompleteness, leading to an under-estimate of the 
early-type population. As a mater of fact, since most of the 
galaxies lost in this mass bin are probably red early-type galax- 
ies, the values measured should be considered as lower limits. In 
order to give a rough estimate of the maximum value, we con- 
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Fig. 9. Stellar mass-color plots in the 2 central Mpc for all the 9 clusters in the sample. Early-type galaxies are shown with red filled 
circles and late-type with blue spiral symbols. Solid line indicate the best fit to the red-sequence and dashed lines 1.5<x limits (see 
text for details). 



sidered that the 50% of galaxies are E/SOs and estimated how 
the fraction of these galaxies would rise at every position (see 
Fig. rTTT >. We observe that, even if we loose half the early-type 
population at all distances, it would remain below the late-type 
population. 

5.3. Morphology-density relation 

Figure[T2]shows the T - £ relation for all the clusters in the sam- 
ple. In all clusters, the fraction of early-type galaxies increases 
monotonically with local surface density, with some heterogene- 



ity due to cluster-to-cluster variations. This relation is not en- 
tirely correlated with the decrease of the early-type fraction with 
radius (the T - R relation) even if there are some similarities. In 
practice, we see that, beyond ~ IMpc, the correlation between 
radius and density (Fig. [8} is weaker with a high dispersion in the 
number density so that we span more different environments. 

Therefore, as shown in previous works, studying the prop- 
erties of galaxies in the center (d < IMpc) and in the out- 
skirts (d > 2Mpc) separately, helps separating different phys- 
ical phenomena. While at the center, the effects of local envi- 
ronment (density) and cluster potential (distance to the center) 
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Table 3. Composition of the red sequence in the 2 central MegaParsecs. Table shows the fraction of galaxies of each morphological 
and mass bins which belong to the red sequence. 



log(M/M Q ) > 11.3 10.3 < log(M/M ) < 11.3 9.5 < log(M/M ) < 10.3 Total 

E/SO 45% 27% 4% 76%~ 

S/I 8% 12% 4% 24% 



are very well correlated, this is not the case in the outskirts. 
Consequently, differentiating the analysis of galaxy properties 
in the outskirts with respect to radial distance R or local number 
density S will tell us if galaxies are entering the cluster in groups 
and evolving afterwards or if they are acquiring the relation be- 
tween morphology and environment/radius once they enter the 
cluster. 

Figure [13] indeed shows the trend with density of the early-type 
fraction averaged over all the clusters for galaxies in the outer 
parts (d > 2Mpc) and in the central regions (d < IMpc). The 
fraction of early-type galaxies is systematically higher in all the 
environments in the center than in the outskirts. This rise of 
the global fraction of E/SO objects reflects the fact that galaxies 
are experiencing a strong morphological segregation in cluster 
centers. In the outskirts of the cluster, the global morphologi- 
cal mixing is correlated with local density, i.e. the fraction of 
early-type galaxies rises with increasing density and the T — E 
relation is preserved. This seems to indicate that morphological 
transformation of galaxies entering the cluster is strongly influ- 
enced by the local environment and more precisely by interac- 
tions with galaxies in their close surroundings. The correlation 
is also present in the cluster center, but with a higher fraction of 
E/SO galaxies at all redshifts. Galaxies therefore continue evolv- 
ing inside the cluster and migrating to the red-sequence but keep- 
ing a memory of their surroundings. 

Are all galaxies evolving at the same speed and independantly 
of their mass? In figure [14] the fractions of each morphologi- 
cal type for three stellar mass ranges in the center and in the 
outskirts of the cluster are shown. We detect different behav- 
iors depending on the considered mass range. Massive galaxies 
log(MJM Q ) > 11.3 are mostly elliptical (~ 70%) in the cluster 
outskirts at all the environments. The fraction is slightly lower 
than the one measured in the cluster center (~ 80%) but remains 
high. This seems to indicate that an important fraction of galax- 
ies in this mass range are already formed when they enter the 
cluster at z ~ 0.5 and have already experienced morphological 
evolution in t heir local environmen t (tidal interactions and merg- 
ing in groups, iBamford et aT]|2009t) . 

Intermediate mass objects (10.3 < log(M*/M Q ) < 11.3) show 
a different trend. If we consider the cluster outskirts, we do 
see a relation between the morphological mixing and the local 
density: galaxies in the higher density bin (log(DlQ) > 2.0) 
are mostly early-type (~ 60%) while the fraction decreases to 
~ 40% at lower densities. This population is therefore experi- 
encing an evolution driven by the local environment when enter- 
ing the cluster and dominates the general trend seen in the whole 
sample (fig. [TBI. If we now look at the central regions, the cor- 
relation is almost absent and the fraction of early-type galaxies 
rises in all the density bins. This again suggests that this popu- 
lation is evolving as it moves to the cluster center but galaxies 
still preserve a memory of their location when they entered the 
cluster. 

If we consider now low massive galaxies (9.5 < log(M t /M ) < 
10.3), we again detect a different behavior. First, this population 



is dominated by late-type objects at almost all densities, even if 
we still measure a higher fraction of E/SO galaxies at high densi- 
ties, i.e. the cross-over between late-type and early-type systems 
occurs at higher density than for the more massive galaxies. It 
seems that morphological transformation is s till going on even 
close to the cluster center dHolden et al.ll2Q07D . Once again, this 
result might be a consequence of incompleteness in this mass 
bin, which leads to an under-estimate of the red early-type pop- 
ulation. As for section I5T21 we show with vertical arrows in fig- 
ure [14] how the fractions would increase if we consider that we 
loose 50% of the early-type population. Taking into account this 
effect reduces the trend, but there are still more late-type galax- 
ies at all densities than than in the more massive bins. 
Therefore, the transition from the blue-cloud to the red sequence 
is performed at different speeds mainly driven by the local en- 
vironment: more massive galaxies drift to the red-sequence at 
earlier epochs. Massive galaxies have converted the gas into 
stars earlier probably because at high redshift their star for- 
mation was not effectively self-limited by feedback. Most of 
these galaxies are already massive, gas-poor and with an early- 
type morphology when they enter the cluster. The transition for 
intermediate/low-mass galaxies is still going-on at z ~ 0.5 and 
seems to be mainly driven by the local environment as stated by 
the correlation between morphology and density at all the clus- 
ter locations. Starbusts resulting from galaxy merging or starva- 
tion processes can be responsible of this transition: galaxy en- 
counters can significantly accelerate the star formation and pro- 
voke the transition to the red-sequence. Semi-analytical models 
show however that this process is effective mainly for midsize 
galaxies which strike the best comprom ise between cross section 
and abundance dMaulbetsch et alj200 7). Galaxy interactions can 
also trigger the AGN phase. The AGN feedback can then con- 
tribute to suppress any residual star formation. 

6. Summary and conclusions 

We have presented a morphological study of a sample of galaxies 
in 9 X-ray luminous clusters at intermediate redshift (z ~ 0.5). 
Wide field imaging provides information up to 5 Mpc from the 
cluster center, well above the virial radius and allows to explore 
the distant outskirts of the clusters, where the local environment 
is dominant with respect to cluster distance. Using photomet- 
ric redshifts we are able to select cluster members with less 
than 20% of field contamination, thanks to a detailed compar- 
ison with spectroscopic redshifts when available. Galaxies are 
then separated into two morphological classes with our auto- 
matic galaxy classifier g alSVM , developped and validated by 
iHuertas-Companv et alj d2008l) . Direct comparison with HST 
imaging of a sub-sample of galaxies in two independant clus- 
ters demonstrates that our method, based on ground-based data, 
is in agreement with HST results for more than 90% of the galax- 
ies up to a magnitude limit r' =23. This magnitude cut which 
defines our working sample is about 1.5 mag. fainter that a L* 
galaxy at z ~ 0.5 and ensures multi-color photometric complete- 
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Fig. 10. Fraction of morphological types (r < 23) as a function of cluster radius. Red points: E/S0 galaxies, blue points: S/l galaxies. 
The total number of objects is shown in each plot. 



ness. In addition, multi-color photometric information, available 
in 4 bands (g', r' , i' and z'), enables an estimate the stellar mass 
of the galaxies from their SED fitting. We estimate that our sam- 
ple is 50% complete down to a stellar mass of 3 10 9 M Q (or 
log(M*/M ) ~ 9.5. All in all, our full galaxy sample contains 
3861 galaxies inside 9 clusters. 
The main results can be summarized as follows: 

- In all clusters, we clearly identify the red sequence of early- 
type galaxies and we are able to build a stellar mass versus 
color sequence of galaxies. Both sequences are mostly pop- 
ulated with massive (log(M,/M ) > 11.3) early-type galax- 



ies. As found by previous studies, the galaxy red sequence 
is already formed at z ~ 0.5 an d early-type gala xies prop- 
erties are already well in place (iPatel et al.ll2009l) . Physical 
processes which occur in dense environments have conse- 
quently a rapid and early influence of the evolution of galaxy 
properties (harassment, ram pressure strippring, etc ...). They 
are very efficient in cluster centers and operate early in the 
cluster histories so that galaxies evolve quickly towards gas 
poor massive early-type galaxies. 
- The population of massive galaxies (log(M,/M ) > 11.3) 
does not seem to have evolved significantly since z ~ 0.5. In 
fact, the fraction of early-type galaxies among this popula- 
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Fig. 11. Morphological mixing as a function of the cluster radius, averaged for all the clusters in the sample and for different mass 
ranges. Filled circles are early-type galaxies and spiral symbols are late-type galaxies. The left panel shows the results for the whole 
sample. Notice that for the low mass bin (top right) the two first points are not significative since there are too few objects. The 
total number of objects is shown in each plot. In the low mass bin, vertical arrows show an estimate of the possible effects of 
incompleteness (see text for details). 




X ■ § L. 

1 : -, 



0.5 1.0 1.5 2.00.0 0.5 1.0 1.5 2.00.0 0.5 1.0 1.5 
Log(D10 [Mpc -2 ]) 



Fig. 12. Morphology-density relation for all the clusters in the 
sample. Blue lines are late-type systems and red lines are early- 
type. The total number of objects is shown in each plot. 



tion is very high and remains almost constant (T-R relation), 
from the cluster center (~ 90%) to the outskirts (~ 70% at the 
virial radius and above), suggesting that these galaxies expe- 
rience no evolution. They were probably formed at higher 
redshifts before entering the cluster in infalling groups and 
their location in the cluster center might be mostly explained 
by a segregation effect. 

When we consider the whole sample, the T - £ relation 
shows a clear correlation between morphology and local 
density with early-type galaxies domi nating the dens est re- 
gions. This was initially pointed out by dDressleifl9 80) in the 
local universe and confirmed later on more distant clusters 



(e.g. lPostman et al1 l2005). The T - S relation is also clearly 
detected in the cluster outskirts, where cluster radius and lo- 
cal density are no longer correlated. It seems therefore that 
the morphological evolution is mainly driven by the local 
environment through galaxy-galaxy interactions (merging, 
starvation), independently of cluster properties. Galaxies en- 
tering the cluster keep a memory of their surroundings and 
evolve accordingly. They evolve however at different speeds, 
depending on the stellar mass. 

Most of the evolution within clusters at z ~ 0.5 is tak- 
ing place in galaxies with intermediate mass (10.3 < 
log(M»/M ) < 11.3) or even lower. These galaxies do show 
a relation between the morphological mixing and the local 
density in both the center and the outskirts with differences 
between dense and loose global environments. In dense envi- 
ronments, early-type galaxies are in place in a similar way as 
for the most massive ones, while in external areas, the T - 2 
relation is similar to the one in the field. So the morpho- 
logical dichotomy is prominent in structures representative 
of galaxy groups but in poorer environments, the morpho- 
logical transformation from spiral-types to early-types has 
not produced significant changes in the galaxy distribution. 
Galaxies in this mass range are experiencing a migration 
from the blue-cloud to the red-sequence as they move into 
the cluster, mainly driven by interactions with their surround- 
ings. The migration is accelerated when galaxies, within 
their local environment, cross the densest parts of the clusters 
in the central Mpc or so. 

Galaxies in the smallest mass bin (9.5 < log(M,/M Q ) < 
10.3) are basically late-type galaxies at all the clusters lo- 
cations and environments (even if they present a trend with 
the local density), suggesting that they do no have started 
the morphological transformation at z ~ 0.5. However, our 
conclusions are weaker for this population because the num- 
ber of galaxies is small, especially in the cluster center and 
results can be affected by incompleteness specially in the 
early-type population. However, even considering that 50% 
of early-type galaxies are lost in this mass bin, this popula- 
tion seems to be dominated by late-type galaxies. 
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Fig. 14. Mean fractions of early-type (filled circles) and late-type (spiral symbols) galaxies as a function of local density in the 
cluster center (top panels) and in the cluster outskirts (bottom panels) and for different mass ranges. Colors indicate three different 
mass ranges: red (log(MJM G ) > 1 1.3), green (10.3 < log(MJM G ) < 1 1.3) and blue (9.5 < log(MJM Q ) < 10.3). The total number 
of objects is shown in each plot. In the low mass bin, vertical arrows show an estimate of the possible effects of incompleteness (see 
text for details). 



The picture which seems to arise from all these results is that 
massive galaxies moved to the red-sequence well before z ~ 0.5. 
Lower mass galaxies (log(M„/M (3 ) < 1 1 .3) are still in process of 
migration from the blue-cloud to the red-sequence. The migra- 
tion is mainly driven by galaxy-galaxy interactions which cause 
a progressive quench of the star formation and/or provoke strong 
starbusts and/or trigger an AGN phase. Only in the central Mpc 
or so, the cluster potential seems to play an important role by 
accelerating the process of migration. 

In order to better trace the general scheme of the morpho- 
logical evolution of galaxies in different environments, two ap- 
proaches are necessary. First, we wish to analyse in more details 
the dependence of these global tendancies with the general prop- 
erties of the clusters themselves. This will be done in a future 
work as soon as the mass distribution of the clusters is available 
using weak lensing mass maps (work in progress, Foex et al. in 
preparation) and/or X-ray properties such as luminosity and tem- 
perature. Clusters of galaxies are not a fully homogeneous popu- 
lation and some of the evolution properties of th eir cluster mem- 
bers a re attached to their intrinsic caracteristics (Po ggianti et al.l 
120091) . Second, it is also necessary to extend this approach to 
more distant clusters to connect clusters at different epochs and 
trace an evolutionary track. In order to better probe the stel- 
lar mass and the morphology traced by the old stellar popula- 



tion, near-infrared imaging is preferred (iHuertas-Companv et al.l 
120091) . 
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